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1 ABSTRACT 

V^Q . Our understanding of how the Galaxy was formed and evolves is severely hampered 

by the lack of precise constraints on basic stellar properties such as distances, masses, 
and ages. Here, we show that solar-like pulsating red giants represent a well-populated 
class of accurate distance indicators, spanning a large age range, which can be used 
to map and date the Galactic disc in the regions probed by observations made by the 
CoRoT0 and Kepler space telescopes. When combined with photometric constraints, 
£NJ ■ the pulsation spectra of such evolved stars not only reveal their radii, and hence dis- 

tances, but also provide well-constrained estimates of their masses, which are reliable 
proxies for the ages of the stars. As a first application we consider red giants observed 
by CoRoT in two different parts of the Milky Way, and determine precise distances 
1 for ~2000 stars spread across nearly 15,000 pc of the Galactic disc, exploring regions 

which are a long way from the solar neighbourhood. We find significant differences in 
the mass distributions of these two samples which, by comparison with predictions of 
synthetic models of the Milky Way, we interpret as mainly due to the vertical gradient 
in the distribution of stellar masses (hence ages) in the disc. In the future, the availabil- 
ity of spectroscopic constraints for this sample of stars will not only improve the age 
determination, but also provide crucial constraints on age- velocity and age-metallicity 
relations at different Galactocentric radii and heights from the plane. 
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1 INTRODUCTION 

Accurate and precise distances to stars are a fundamental 
cornerstone of astronomy. Such data play a key role in shap- 
ing theories of the history and fate of our Galaxy. Another 
fundamental cornerstone is provided by stellar ages. Ages 
can be estimated only for a limited sample of individual 
field stars, and then only by using either model-dependent 
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techniqu es or empirical m ethods that still need careful cali- 
bration (|Soderblomll2010l h These two cornerstones are es- 
sential to Galactic Archaeology, the study of the chemi- 
cal compositions and dynamic motions of stars of different 
ages, within which is coded information on the origin and 
subsequent evolution of the Milky Way l|Turon et al.ll2008l ; 
iFreeman fc Bland-Hawthornll2002h . There is a modern con- 
sensus among astronomers that a variety of processes play 
a role in shaping our Galaxy, for example gas accretion and 
in-situ star formation, mergers, and dynamical processes in 
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the Galaxy (such as diffusion and migration of stars). Galac- 
tic Archaeology research has now reached a point where 
so-called chemo-dynamical models are needed to quantify 
the relative importance of each of these processes, over the 
course of more than 10 Gyr of evolution of the Milky- Way. 
Stringent observational constraints are in turn desperately 
needed to test these model predictions. 

Particularly important observational constraints are 
the age- velocity and age-metal licity relations at different po- 
sition s i n the Galaxy (see e.g. [Freeman fe Bland-Hawthornl 
120021 ; iMinchev. Chiappini fc Martial l2012f ). These 
relations are uncertain, or limited to neighbour- 
ing stars (e.g. from the Geneva-Copenhagen Survey 
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( Hol mberg. Nordstrom fc Andersenll2009l ; ICasagrande et all 
l201ll ). which is confined to distances less than about 100 
pc from us). Future planned surveys hope to fill this gap 
by combining precise p roper motions from the ESA-GAIA 
ijPerrvman et al.l l200lT l astrometric mission with massive 
spectroscopic follow-up from ground-based telescopes, to 
extend the region for which distances and chemical compo- 
sitions are known out to as far as 10,000 pc. However, this 
will not happen before the beginning of the next decade, 
and it will still need to confront the large and systematic 
uncertainties affecting the age determinations of single 
stars. 

Red giants are cool highly luminous stars, which, by 
virtue of covering a wide domain in mass, age, chemical 
composition, and evolutionary state, are in principle an im- 
portant source of data for testing chemo-dynamical models. 
Previously, age estimation of red giants has had to rely on 
constraints on their surface properties only. A severe lim- 
itation of such an approach is that during the red-giant 
phase stars of significantly different age and distance end 
up sharing very similar observed surface properties making 
it extremely hard to discriminate evolutionary states of field 
giants belonging to the Galactic-disc population. That situ- 
ation has now changed thanks to asteroseismology, with the 
detection of solar - like pulsations in thousand s of red giants 
foe Ridder et al.1 12009] ; iBedding etHI l20ld; iHekker et all 
120111 ) observ ed by the C0R0T (iBaglin fc Fridlundl 120061 ) 
and Kepler jGilliland et al.l 120101 7 space telescopes. The 
pulsation frequencies may be used to place tight con- 
straints on the fundamental stellar pr operties, inc l uding 
radius, mass and e v olutionary state ( Stello et alj 2008 ; 
Kallinger et al.1 12010|; IMosser et~ai1 l20ld ; iMontalban et al.1 
2010l ; [Bedding et al.ll201ll), the properties of Helium ionisa- 



tio n regions llMiglio et alj|20ld ). and internal rotation (see 
e.g. lBeck et al.ll2012l ). Here, we show it is possible to use pul- 
sating red giants as accurate distance indicators, and tracers 
of stellar population ages. 



2 METHOD 

Radii and masses of solar-like oscillating stars can be es- 
timated from the average seismic parameters that charac- 
terise their oscillation spectra: the so-called average large 
frequency separation (A//), and the frequency correspond- 
ing to the maximum observed oscillation power (^ ma x). 

The large frequency spacing is predicted by theory to 
scal e as the square root of the mean density of the star (see 
e.g. IVandakurovIl 19671 ; lTassoul|[l980h : 



M/Mq 



(1) 



where Auq = 135 /iHz. The frequency of maximum power is 
expected to be proportional to the acoustic cutoff frequency 



expected to b e proportio nal to the acoustic cutoft frequency 
llBrown et al.lll99ll; iKieldsen fc Bedding||l995l ;l Mosser et al.l 
|2010| ; iBelkacem et al.ll201ll ). and therefore: 
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where ^ max ,0 = 3100 /iRz and T eff , = 5777 K. 

When no information on distance/luminosity is avail- 
able, which is the case for the vast majority of field stars ob- 
served by C0R0T and Ke pler, Eq. [T]and[2]may be solved to 
derive M and R (see e.i dStello et aT] |2008; Kall inger et al.l 



derive M and H (see e.g . 
l20ld : IMosser et all boiO) 
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We determined stellar radii and masses by com- 
bining the available seismic parameters f max and Av 
with effective temperatur es r e g. The l a tter were de- 
termined using 2MASS |Skrutskie et al.1 l2006h J and 
K 3 photometry and the colou r - T e ff calibrations by 
lAlonso. Arribas fc Martmez-Rogerl (|l999l ). which depend 
only weakly on metallicity. 2MASS colours were trans- 
formed into the CIT photome t ric s ystem used by 
lAlonso. Arribas fc Martmez-Rogerl (Il999 | ) using the rela- 



tion s available in Alonso. Arribas fc Martinez-Rogerl (Il998l ) 
and lCarpenterl <|200ll ). 

We then computed luminosities L using the Stefan- 
Boltzmann law L — 4,-KR 2 aT* B , where a is the Stefan- 
Boltzmann constant. The distance modulus of each star was 
determined as K' s0 — K 3 o, where K a o is the de-reddened ap- 
parent 2MASS K s magnitude, and K' a0 the absolute K s 
magnitude. The latter was obt ained by combining L and 
the bolometric corrections from lGirardi et all l|2005h . which 
are based on the Castelli & Kurucz (2005) ATLAS9 model 
atmospheres, for the range of T c s and log g under consider- 
ation. 

We took into account the effect of interstellar extinc- 
tion on the magnitude and colour of each star using the 
3D model of Galactic extinction in the V b a nd (^4 y) by 
iDrimmel. Cabrera-Lavers fc Lopez-Corredoiral |2003D- The 
extinctio ns in the J and K s band passes were determined 
following iFiorucci fc Munaril (|2003| j assuming the spectral 
energy distribution of a Kl giant. Since the extinction is 
distance dependent, we iterated the procedure until the de- 
rived distance does not vary by more than 1%. We chose to 
consider magnitudes in the near-IR to reduce the effect of 
interstellar reddening in both the determination of T e g and 
apparent de-reddened magnitudes. 

By analogy with the well-known period-luminosity re- 
lation u sed to estimate distances of classical pulsators 
fsee e.g. iFeast fc Walker! Il987l : iMadore fc Freedmanl Il99ll : 
iBono et alj|2010l . for a review), we can explicitly write a re- 
lation between the distance and the pulsation properties of 
solar-like oscillators. In contrast to the case of single-mode, 
radial pulsators we can estimate stellar radii without making 
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any assumption^] on mass-metallicity-luminosity relations. 
Hence we find: 



logd= l + 2.51of 
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where d is expressed in parsec, m^ \ is the apparent bolo- 
metric magnitude of the star, and Mbol,© the absolute solar 
bolometric magnitude. Although this expression is model- 
independent, it builds upon Eq[T]and Eq.[2]which are based 
on simplifying assumptions that need to be independently 
verified, and supported with empirical tests. 

Tests against independent estimates of mass and ra- 
dius suggest that R and M derived using average asteroseis- 
mic param eters are accurate to within 5% and f0%, respec- 
tively (see iMiglio et al. 20f 2 for a review on the te s ts per - 
formed so fa r seelMiglid (|2012l l : ISilva Aguirre et all <|2012T) ; 
iHuber et al.l (|20ia T for tests on nearby stars). 



the apparent de-reddened luminosity I. Since d oc (L/l) 1 ^ 2 , 
the overall effect of reddening on the distance itself is partly 
reduced. By performing 500 realisations of the data assum- 
ing the uncertainties described above to be Gaussian, we 
find a median intrinsic uncertainty of 5% in distance, 3.6 % 
in radius, and 10 % in mass. To explore th e effect of possible 
syste matic offsets in the T e ff scale (see e.g. ICasagrande et all 
2010), we increased/decreased T e g by 100 K, and found that 
the distance estimate was affected by 2.5% (the radius by 
1% and the mass by 3%). This leads to an overall uncer- 
tainty on the distance of 5% (random) + 5% (systematic 
due to radius determination via scaling relations) + 2.5% 
(systematic due to the T c ff scale). 

The positions of the giants in the LRcOl and LRaOl are 
illustrated in Fig. [T] 



3 DATA AND UNCERTAINTIES 

The CoRoT photometric time series used in this work 
were obtained in the so-called exofield during the first 
long CoRoT runs: LRcOl (I = 37°, b = -7°) and LRaOl 
(I = 217°, b = —2°). These long runs lasted approx- 
imately 140 days, providing us with a frequency resolu- 
tion of a bout 0.08 ^tHz. A first a nalysis of these data was 
made by |Pe Ridder et all l|2009l ). iHekker et all (|2009l ) and 
iKallinger et all l|2010l ). Solar-like oscillations were detected 
in 435 and 1626 giants belonging to LRaOl and LRcOl, 
respectively. To derive M,R, and distance we use the val- 
ues (and unc e rtaint ies) of i/ m „ and Au as determined by 
iMosser et al] |20ld ). The typical uncertainty on tWx and 
Au for the 150-d long CoRoT observations is of the order of 
2.4% and 0.6%, respectively. 

The uncertainties on the stellar properties are estimated 
using Monte Carlo simulations adopting the following con- 
straints: 

• Apparent magnitudes: uncertainties on J and K a were 
taken from the 2MASS photometry available in the EXO- 
DAT (|Deleuil et al|2 009) catalogue (the median uncertainty 
is 0.02 mag). 

• Extinction/reddening: a random error in Ay of 0.3 mag 
was considered both in determining T c g photometrically and 
in de-reddening K a apparent magnitudes. 

• T c ff. for each star we considered two sources 
of uncertainty. The first (100 K) was due to un- 
certainties on the colour-T e ff ca libration itself 
l|Alonso. Arribas Sz Martinez- Roger! 119991 ). The second 
one was due to uncertainties on reddening. This resulted 
in a median combined uncertainty on T e g of ~ 190 K 
(calibration+reddening) . 

When determining the distance, reddening affects not 
only the de-reddened apparent bolometric magnitude, but 
also the determination of L (mostly via T e ff)- A higher Ay 
increases the estimated T e fj, hence L, but it also increases 



1 Excluding the weak dependence of bolometric correction on 
metallicity. 



4 DIFFERENTIAL POPULATION STUDY 

Before comparing populations of giants observed in the two 
fields, we checked whether different biases were introduced in 
the target selection in LRaOl a nd LRcOl. We retrieved from 
EXODAT (|Deleuil et all 120091 ') photometric information on 
all the stars in the field, as well as the targets observed. Tar- 
gets within each field of view were selected largely on the 
basis of colour-magnitude criteria. Within the observed tar- 
gets, solar-like oscillations were searched for in stars belong- 
ing to a limited col our-magnitude dom ain: 0.6 < J — K s < 1 
and K s < 12 fsee IMosser et al.ll201(t) . Restricting to this 
domain, we find no significant difference in the target se- 
lection bias applie d to LRcOl compared to LRaOl (see also 
IMiglio et aHl2012r ). 

We then compare the distributions of radius and mass 
of stars in LRaOl and LRcOl. Since scaling relations were 
not tested at high luminosities we excluded from the sample 
the few stars with L > 200 L©, and stars with no 2MASS 
photometry available. The fraction of targets excluded rep- 
resents < 3% of the whole sample. The distributions of R 
and M are reported in Fig. [2] 

We applied the Kolmogorov-Smirnov (K-S) test to the 
distributions of mass and radius to quantify differences be- 
tween the two populations. When comparing radii, the null 
hypothesis (LRaOl and LRcOl samples are drawn from the 
same parent distribution) cannot be rejected, or the differ- 
ence between the two populations is marginally significant 
at best. On the other hand, we find the difference between 
the mass distribution of the two populations to be highly 
significant (K-S probability higher than 99.9 %). 

To test the impact on the K-S test of statistical uncer- 
tainties on mass and radius, we perturbed the observed radii 
and masses adding random offsets drawn from a Gaussian 
distribution having a standard deviation equal to the esti- 
mated uncertainty on the mass / radius. We generated 1000 
realisations and performed a K-S test on each realisation. 
The distribution of results from the 1000 K-S tests confirms 
that while comparing radius the null hypothesis cannot be 
rejected, the difference in the mass distribution is significant 
(in 95% of the realisations the K-S probabilities are higher 
than 99%) . 
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Figure 1. Positions of red giants in the two observed fields (CoRoT designations LRcOl and LRaOl) projected on the Galactic plane 
(right panel) and on a plane perpendicular to the disc (left panel). 
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Figure 2. Radius and mass distributions of stars in LRaOl (blue) 
and LRcOl (red) obtained using the method described in Sec. [2] 



4.1 Comparison with synthetic stellar populations 

For a meaningful comparison between observed and syn- 
thetic populations, we applied to the synthetic population 
the same selection criteria based on colour and magnitude 
adopted in the target s electi on. Moreover, as discussed in 
detail in iMosser et alj (HpiQ), since no significant bias in 
the distribution of the targets is present in the ^ ma x range 
between 6 fiRz and 80 jiRz, we only considered observed 
and simulated stars with ^ ma x in this frequency range. We 
considered synthetic populations compu ted with two codes : 
the Besancon Model o f the Milky Way jRobin et al.ll2003 ). 
and with TRILEGAL jGirardi et al.ll2005t ). 

Simulations with both TRILEGAL (see Fig. [3} and the 
Besancon Model show that, although a similar distribution 
of radius is expected in LRaOl and LRcOl, the age (hence 
mass) distribution of the two populations of giants is ex- 



pected to differ. Stars in LRcOl are expected to be older, on 
average, than those in LRaOl. Additional tests presented in 
iMiglio et al.1 (|2012l ) show that, in the simulations, the dif- 
ferent mass distribution expected in the two fields is indeed 
related to the assumed age-dependent vertical scale height 
of the thin disc. Our aim here is not to find the synthetic 
population that best matches the data, but to show that the 
difference we see in the observed distribution is in qualitative 
agreement with the simulations which include an increase of 
the disc scale height with age. 

The estimated masses of our sample provide important 
constraints on the stellar ages. Once a star has evolved to the 
red-giant phase, its age is determined to good approxima- 
tion by the time spent in the core-hydrogen burning phase, 
and this is predominantly a function of mass. The CoRoT 
giants cover a mass range from ~ 0.9 to ~ 3 M©, which 
in turn implies an age range spanning ~ 0.3 to ~ 12 Gyr, 
i.e., the entire Galactic history. To estimate the age of gi- 
ants in the two observed populations we use PARAM, a 
Bayesian stellar parame t er est imation method described in 
detail i n l da Silva et al] (120061 ). For this work we adapted 
paramEI to include as additional observational constraints 
i/ mal , Au and, when available, the evolutionary status of the 
star (core-Helium or hydrogen-shell burning phase), as de- 
ter mined by the period spacing of gravity-dominated modes 
(see lMosser et al1l201ll ). Having knowledge of the evolution- 
ary status is particularly useful in the age determination of 
stars with an estimated radius typical of RC giants. Since 
the mass of core-He burning stars is likely to be affected by 
mass loss, the age-mass relation in red-clump stars may be 
differ ent from that of RGB stars (see e.g. Fig. 1 in iMiglid 
120121 ). 

We applied this parameter estimation method to the 
sample of giants observed in LRcOl and LRaOl and obtained 
distributions of age, mass and radius. The estimated uncer- 
tainty on the age is of the order of 30-40% given that no 
information on the metallicity is available, hence a broad 
prior on [Fe/H] is assumed . The results obtained with a 
model-based approach support the interpretation that the 



2 This extended version of PARAM will be made 
available, via an interactive web form, at the URL 
http : //stev. oapd. inaf . it/param 
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Figure 3. TRILEGAL simulations. Distributions of radius (left panel), mass (middle panel), age (right panel) of giants in synthetic 
populations computed with TRILEGAL for stars in the field LRaOl (blue) and LRcOl (red). Stars in the synthetic population were 
selected in magnitude (K s < 12), colour (0.6 < J — K s < 1), and ^ max (6 /^Hz < l>max < 80 fiHz) to account for target selection effects. 



differences between the observed population in LRaOl and 
LRcOl are due to a different mass, hence age, distribution. 



5 CONCLUSIONS 

In this paper we have shown that solar-like oscillating gi- 
ants are key tracers of stellar populations in the Milky Way. 
When combined with photometric constraints, the pulsa- 
tion spectra of solar-like oscillating giant stars not only re- 
veal their radii, and hence distances, but also provide well- 
constrained estimates of their masses, which are reliable 
proxies for the ages of the stars. 

We considered red giants observed by the CoRoT space 
telescope in two regions located at different positions on the 
sky, with CoRoT designations LRcOl (I = 37°, b = -7°), 
and LRaOl (I = 217°, b = -2°). We found a significant dif- 
ference in the mass distributions of the populations of the 
stars in the two fields, with stars in LRcOl having a lower 
average mass than those in LRaOl. Differences in radius are 
in contrast marginal at best. To interpret these findings we 
have compared the observed distributions with the distribu- 
tions given by synthetic population simulations of red giants 
for the same fields of the Galaxy covered by the observations. 
The synthetic populations show a difference in the mass dis- 
tributions that agrees qualitatively with that found in the 
observed populations. 

On the basis of these comparisons we interpret the 
differences in the mass distributions as being due mainly 
to the different average heights of the observed fields be- 
low the galactic plane. Since it is believed that dynam- 
ical processes in the disc increase the velocity dispersion 
of stars with time, it follows that older stellar populations 
reach greater heights above and below the plane. What we 
observe with CoRoT, i.e. that the field higher below the 
plane has a larger fraction of old (i.e low mass) stars, thus 
seems to support this theoretical expectation. Other data 
may hint at a similar dependence, as for instance an in- 
creas e of the velocity dispersion of stars wi th increasing 
age l|Holmberg. Nordstrom fc Andersen] 120091). or the cor - 
relation between scale height and [O /Fe] (|Bovv et al.ll2012l ) , 
but they are only available for stars in the solar vicinity. Fu- 
ture analyses of pulsating red giants sampling the Milky Way 
at different radii and heights from the plane, when comple- 
mented with their chemical abundances, will set even tighter 



constraints and enable to fully disentangle both the radial 
and vertical structure of the Galactic disc. 

The findings presented here provide an example of the 
detailed picture of Galactic populations inferred from pul- 
sating red giants. We can map regions that are more than 
a factor ten further away in distance compared to what 
has previously been possible, and we can explore (with a 
homogeneous sample) a wide age interval sampling look- 
back times as long as the age of the Galaxy. This method, 
when applied to the various regions explored by CoRoT and 
Kepler, and complemented by chemo-dynamical constraints 
from spectroscopic analyses, will provide the gold standard 
for current and future surveys of the Milky Way. 
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